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Abstract. A new scheme for testing the nuclear matter 
(NM) equation of state (EoS) at high densities using con- 
straints from compact star (CS) phenomenology is applied 
to neutron stars with a core of deconfined quark matter 
(QM). An acceptable EoS shall not to be in conflict with 
the mass measurement of 2.1 ± 0.2 M Q (1 a level) for 
PSR J0751+1807 and the mass radius relation deduced 
from the thermal emission of RX J1856-3754. Further con- 
straints for the state of matter in CS interiors come from 
temperature-age data for young, nearby objects. The CS 
cooling theory shall agree not only with these data, but 
also with the mass distribution inferred via population 
synthesis models as well as with LogN-LogS data. The 
scheme is applied to a set of hybrid EsoS with a phase 
transition to stiff, color superconducting QM which fulfills 
all above constraints and is constrained otherwise from 
NM saturation properties and flow data of heavy-ion colli- 
sions. We extrapolate our description to low temperatures 
and draw conclusions for the QCD phase diagram to be 
explored in heavy-ion collision experiments. 



1. Introduction 

Recently, new observational limits for the mass and 
the mass-radius relationship of CSs have been obtained 
which provide stringent constraints on the equation of 
state of strongly interacting matter at high densities, see 



Klahn et al. (2006) and references therein. In this latter 



work several modern nuclear EsoS have been tested re- 
garding their compatibility with phenomenology. It turned 
out that none of these nuclear EsoS meets all constraints 
whereas every constraint could have been fulfilled by some 
EsoS. As we will point out in this contribution, a phase 
transition to quark matter in the interior of CSs might 
resolve this problem. In the following we will apply an 
exemplary EoS for NM obtained from the ab-initio rela- 
tivistic Dirac-Brueckncr-Hartree-Fock (DBHF) approach 
using the Bonn A potential ([van Dalen et al. 2005|) . There 
is not yet an ab-initio approach to the high-density EoS 
formulated in quark and gluon degrees of freedom, since 



it would require an essentially nonperturbativc treatment 
of QCD at finite chemical potentials. For some promis- 
ing steps in the direction of a unified QM-NM description 
on the quark level, we refer to the nonrelativistic poten- 



tial model approach by Ropke et al. (1986) and the NJL 
model one by Lawley et al. (2006) Simulations of QCD on 
the Lattice meet serious problems in the low-temperature 
finite-density domain of the QCD phase diagram relevant 
for CS studies. However, there are modern effective ap- 
proaches to high-density QM which, albeit still simplified, 
focus on specific nonperturbativc aspects of QCD. They 
differ from the traditional bag model approach and al- 
low for CS configurations with sufficiently large masses, 
For our QM description we em- 



Alford et al. (2006) 



ploy a three-flavor chiral quark model of the NJL type 
with selfconsistent mean fields in the scalar meson (cou- 
pling Gs) and scalar diquark (coupling Go = fjD Gs) 
channels (Blasch ke et al. 2005[) . generalized by including 
a vector meson mean field (coupling Gy = rjv Gs), see 



Klahn et al. (2006a) 



We show that the presence of a QM core in the in- 
terior of CSs does not contradict any of the discussed 
constraints. Moreover, CSs with a QM interior would be 
assigned to the fast coolers in the CS temperature-age di- 
agram. Another interesting outcome of our investigations 
is the prediction of a small latent heat for the deconfine- 
ment phase transition in both, symmetric and asymmetric 
NM. Such a behavior leads to hybrid stars that "masquer- 
ade" as neutron stars and has been discussed earlier by 



Afford et al. (2005)| for a different EoS. This finding is of 
relevance for future heavy-ion collision programs at FAIR 
Darmstadt. 

2. The flow constraint from HICs 

The behaviour of elliptic flow in heavy-ion collisions is 
related to the EoS of isospin symmetric matter. The 
upper and lower limits for the stiffness deduced from 
such analyses (jDaniclcwi cz et al. 2002[) are indicated in 
Fig. 1 as a shaded region. The nuclear DBHF EoS is soft 
at moderate densities with a compressibility K — 230 
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Fig. 1. Constraint on the high-density behavior of the 
EoS from simulations of flow data from heavy-ion collision 
experiments (shaded area from Daniele wicz et al. 2002[) 
compared to the nuclear matter and hybrid EsoS discussed 
in the text. 



MeV (Ivan Dalen et al. 20041 |Gross-Boelting et al. 1999| , 
but tends to violate the flow constraint for densities above 
2-3 times nuclear saturation. As a possible solution to this 
problem we adopt a phase transition to QM with an EoS 
fixed to sketch the upper boundary of the flow constraint. 
In order to obtain an EoS as stiff as possible we use a 
vector coupling of rjv = 0.50 and a diquark coupling of 
i]£> = 1.03. Herewith the EoS is completely fixed. 

3. Constraints from astrophysics 

3.1. Maximum mass and mass-radius constraints 

These most severe constraints come in particular from the 
mass measurement for PSR J0751+I807 (|Nice et al. 2005[) 
giving a lower limit for the maximum mass w 1.9 Mq at 
1<T level, and from the thermal emission of RX J1856- 
3754 ( jTriimper et al. 2004P providing a lower limit in the 
mass-radius plane with minimal radii R > 12 km. These 
constraints can only be fulfilled by a rather stiff EoS. The 
most stiff quark matter contribution to the EoS which 
still fulfills the flow constraint in symmetric matter cor- 
responds to rjv = 0.5 with a maximum mass for hybrid 
stars W 2.1 Mq, rather independent of the choice of r\r> 
which fixes the critical mass for the onset of deconfine- 
ment, see Figs. 2, 3. For a more detailed discussion, see 
Klahn et al. (2006)[ |Kl"ahn et al. (2006a) 



3.2. Cooling constraints 

Direct Urea (DU) processes are flavor-changing pro- 
cesses with the prototype being n — > p + e~ + v e 
( jGamow and Schoenberg 1941| ) , providing the most effec- 



Fig. 2. Stable CS configurations for neutron stars 
(DBHF) and hybrid stars, characterized by the param- 
eters r\D and rjv of the quark matter EoS. 



2.5 
2.0 
1.5 

S i.o 

0.5 



o, 




DBHF 




T] D = 0.92,T1 V 


= 0.0 


n D = i-oo, ti v 


= 0.0 


n D = i.oo, ti v 


= 0.5 


U D = 1-02, T) v 


= 0.5 


n D = 103, ri v 


= 0.5 



10 



R[km] 



12 



14 



Fig. 3. Mass-radius relations for CSs with possi- 
ble phase transition to deconfined quark matter, see 



Klahn et al. (2006a) 



tive cooling mechanism in the hadronic layer of compact 
stars. It acts if the proton fraction x exceeds the DU 
threshold xdu, % = n p /(n n + n p ) > xdu- The threshold 
is given by xdu = 0.11 (jLattimer et al. 1991|) and rises 
up to xdu = 0.14 upon inclusion of muons. Although 
the onset of the DU process entails a sensible depen- 
dence of cooling curves on the star masses, hadronic cool- 
ing with realistic pairing gaps is not sufficient to explain 
young, nearby X-ray dim objects, like Vela, with typical 
CS masses, not exceeding 1.5 M Q (IBlaschke et al. 2 004, 



Grigorian et al. 2005). The point on the stability curve 
in Fig. 2 marks the DU threshold density for the 
DBHF EoS. Quark matter DU processes provide en- 
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Fig. 4. Cooling evolution for hybrid stars of different 
masses given in units of Mq. Note that Vela is described 
with a typical CS mass not exceeding 1.45 Mq. From 




Popov et al. 2006a 
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Fig. 5. Comparison of observational data for the LogN- 
LogS distribution with results from population synthesis 
using hybrid star cooling according to Popov et al. 200 6a. 



hanccd cooling, characterized by the diquark pairing gaps 
(IBlaschke et al. 20001 |Page et al. 2000[) and their densit y 



dependence ( Grigorian et al. 2005 Popov et al. 2006a). 



For a recent review, see Sedrakian (2007) 



To verify this rather heuristic approach we apply ex- 
plicit calculations of the cooling of hybrid configurations 
which shall describe present data of the temperature- age 
distribution of CSs. The main processes in nuclear matter 
that we accounted for are the direct Urea, the medium 
modified Urea and the pair breaking and formation pro- 
cesses. Furthermore we accounted for the ISo neutron 
and proton gaps and the suppression of the 3P2 neutron 
gap. For the calculation of the cooling of the quark core 
we incorporated the most efficient processes, namely the 
quark modified Urea process, the quark bremsstrahlung, 
the electron bremsstrahlung and the massive gluon-photon 
decay. In the 2-flavor superconducting phase one color of 
quarks remains unpaired. Here we assume a small resid- 
ual pairing (Ax) of the hitherto unpaired quarks. For de- 
tailed discussions of cooling calculations and the required 
ingredients see IBlaschke et al. 20041 Popov et al. 2006a] Q] 
and references therein. The resulting temperature-age re- 
lations for the introduced hybrid EoS are shown in Fig. 4. 
The critical density for the transition from nuclear to 
quark matter has been set to a corresponding CS mass of 
-Merit = 1-22 Mq. All cooling data points are covered and 
correspond to CS configurations with reasonable masses. 
In this picture slow coolers correspond to light, pure neu- 
tron stars (M < M cr ;t), whereas fast coolers are rather 
massive CSs (M > M cr ;t) with a QM core. 

Another constraint on the temperature-age relation is 
given by the maximum brightness of CSs, as discussed by 



observational efforts one has not observed very hot NSs 
(logT > 6.3 - 6.4 K) with ages of 10 3 - 10 4 5 years. Since 
it would be very easy to find them - if they exist in the 
galaxy - one has to conclude that at least their fraction is 
very small. Therefore a realistic model should not predict 
CSs with typical masses at temperatures noticeable higher 
than the observed ones. The region of avoidance is the 
hatched trapezoidal region in Fig. 4. 

The final CS cooling constraint in our scheme is 
given by the Log N-Log S distribution, where N is 
the number of sources with observed fluxes larger than 
S. This integral distribution grows towards lower fluxes 
and is inferred, e.g., from the ROSAT all-sky survey 
flNeuhauser and Triimper 1999D . The observed Log N-Log 
S distribution is compared with the ones calculated in the 
framework of a population synthesis approach in Fig. 5. A 
detailed discussion of merits and drawbacks can be found 



Popov et al. (2006) 



Altogether, the hybrid star cooling behavior obtained 
for our EoS fits all of the sketched constraints under the 
assumption of the existence of a 2SC phase with X-gaps. 

4. Outlook: The QCD phase-diagram 

Within the previous sections we exemplified how to ap- 



ply the testing scheme introduced in Klahn et al. (2006) 



Grigorian (2006) It is based on the fact that despite many 



to the modeling of a reliable hybrid EoS with a NM-QM 
phase transition that fulfills a wide range of constraints 
from HICs and astrophysics. In a next step we extend 
the description to finite temperatures focusing on the be- 
haviour at the transition line. For this purpose we ap- 
ply a rclativistic mean-field model with density-dependent 
masses and couplings ( |Typel 2005D adapted such as to 
mimick the DBHF-EoS and generalize to finite temper- 
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Fig. 6. Phase diagram for isospin symmetry using the 
most favorable hybrid EoS of the present study. The NM- 
2SC phase transition is almost a crossover. The model 
DD-F4 is used as a finite-temperature extension of DBHF. 
For the parameter set (t]d = 0.75, rjv = 0.0) the flow con- 
straint is fulfilled but no stable hybrid stars are obtained. 



atures (DD-F4) . Fig. 6 shows the resulting phase diagram 
including the transition from nuclear to quark matter 
[i)d = 1-030, rjv = 0.50) which exhibits almost a crossover 
transition with a negligibly small coexistence region and a 
tiny density jump. At temperatures beyond T ~ 45 McV 
our NM description is not reliable any more since contri- 
butions from mesons, hyperons and nuclear resonances are 
missing. This will be amended in future studies. 

5. Conclusions 

We have presented a new scheme for testing nuclear mat- 
ter equations of state at supernuclear densities using con- 
straints from neutron star and HIC phenomenology Mod- 
ern constraints from the mass and mass-radius-relation 
measurements require stiff EsoS at high densities, whereas 
flow data from heavy-ion collisions seem to disfavor too 
stiff behavior of the EoS. As a compromise we have pre- 
sented a hybrid EoS with a phase transition to color su- 
perconducting quark matter which, due to a vector meson 
meanfield, is stiff enough at high densities to allow com- 
pact stars with a mass of 2 M . Such a hybrid EoS could 
be superior to a purely hadronic one as it allows a faster 
cooling of objects within the typical CS mass region. This 
way, young nearby X-ray dim objects such as Vela could be 
explained with masses not exceeding 1.5 Mq. The present 
hybrid EoS predicts hybrid stars that "masquerade" as 
neutron stars, suggesting only a tiny density jump at the 
phase transition. This characteristics is also present for 
the symmetric matter case and persists at higher temper- 



atures in the QCD phase diagram. It is suggested that the 
CBM experiment at FAIR might softly enter the quark 
matter domain without extraordinary hydrodynamical ef- 
fects from the deconfinemcnt transition. 
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